Fermi and Non-Fermi Liquid Behavior of Quantum Impurity Models:* 
A Diagrammatic Pseudo-particle Approach 
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We review a systematic many-body method capable of describing Fermi liquid and Non-Fermi 
liquid behavior of quantum impurity models at low temperatures on the same footing. The crossover 
to the high temperature local moment regime is covered as well. The approach makes use of a pseudo- 
particle representation of the impurity Hilbert space in the limit of infinite Coulomb repulsion U as 
well as for finite U . Approximations are derived from a generating Luttinger-Ward functional, in 
terms of renormalized perturbation theory in the hybridization V . Within a "conserving T-matrix 
approximation" (CTMA), including all two-particle vertex functions, an infinite series of leading 
infrared singular skeleton diagrams is included. The local constraint is strictly enforced. Applied to 
the SU(N) X SU(M) multichannel Anderson model the method allows to recover the Fermi liquid 
behavior of the single channel case, as well as the non- Fermi liquid behavior in the multi-channel case. 
The results are compared with the "non-crossing approximation" (NCA) and with data obtained 
by the numerical renormalization group and the Bethe ansatz. The generalization of the method to 
the case of finite on-site repulsion U is presented in a systematical way and solved on the level of a 
generalized NCA, fully symmetric with respect to all virtual excitations of the model. 



I. INTRODUCTION 



Over the past two decades the problem of correlated 
electrons on a lattice has emerged as a central theme 
of condensed matter theory. Since the mid seventies 
numerous metallic systems have been found in which 
the conduction electron system is believed to be domi- 
nated by a strong on-site Coulomb iiitjeraction. Exam- 
ples are the heaw fermion compoundstl and the cuprate 
superconductorsH, but systems like the two-dimensional 
electron gas in quantum Hall deviceaj are governed by 
similar physics. Typical model Hamiltonians are the 
Hubbard model and the t-J model for the single-band sit- 
uation, and the Anderson or Kondo lattice models in the 
multi-band case. In the event that the local Coulomb re- 
pulsion in these models exceeds the band width, conven- 
tional many-body theory, i.e. perturbation theory (in- 
cluding infinite resummations) in the Coulomb repulsion 
does not work anymore. The obvious alternative, per- 
turbation expansion in terms of the kinetic energy meets 
with severe difficulty, mainly due to the infinite degener- 
acy of the ground state of a lattice model in the limit of 
zero hopping (a collection of independent atoms) and the 
non-canonical commutation relations of the field opera- 
tors of the eigenstates in the atomic limit. It is therefore 
desirable to develop new methods which are capable of 
dealing with these problems. 

To make it clear what is the principal obstacle for ap- 
plying the well-established many-body techniques to this 
class of problems let us consider the Hubbard model with 
large on-site repulsion [/, where each lattice site can ei- 
ther be empty (state | 0)), singly occupied (|t), |J,) for 
an electron with spin projection 1,1) or doubly occu- 



pied (I 2)). The dynamics of an electron will be very 
different according to whether it resides on a singly or 
doubly occupied site. For less than half band filling and 
large repulsive U the doubly occupied states will be very 
costly in energy, and will contribute to the low energy 
physics only via virtual processes. This is the typical 
feature of strongly correlated electrons: the dynamics is 
constrained to a subspace of the total Hilbert space. The 
necessary projection onto this subspace is difficult to per- 
form within conventional many-body theory and requires 
the development of new theoretical tools. This article 
deals with the method of auxiliary particles introduced 
to effect the projection in Hilbert space, while keeping 
most of the desirable features of renormalized perturba- 
tion theory. 

The idea of auxiliary particles is simple: for each of 
the four states | 0),|t),|i),| 2) at a given lattice site 
(considering one orbital per site), a particle is introduced, 
which is created out of a vacuum state without any lattice 
site at all. The Fermi character of the electron requires 
that two of the auxiliary particles are fermions, e.g. the 
ones representing \]) and ||), and the remaining two are 
bosons. To ensure that a lattice site is in one of the four 
physical states, the number of auxiliary particles is con- 
strained to be exactly one. Compared to alternative ways 
of effecting the projection in Hilbert space, the auxiliary 
particle method has the advantage of allowing one to use 
the machinery of quantum field theory, provided the con- 
straint can be dealt with in a satisfactory way. As will be 
shown below, the constraint is closely related to a ficti- 
tious gauge symmetry in the Fock space of pseudo parti- 
cles. The formulation in terms of renormalized perturba- 
tion theory in the kinetic energy presented in this paper 
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allows to conserve the gauge symmetry in any chosen con- 
serving approximation. Conserving approximations, on 
the other hand, are generated by functional derivation of 
a Luttinger-Ward functional. The Luttinger-Ward func- 
tional is chosen taking into account information on im- 
portant physical processes as well as using expansion in a 
small parameter, e.g. the hybridization energy in units of 
the conduction electron band width. Auxiliary particle 
representations were pioneered by AbrikosovQ, who first 
introduced a feunionic representation for local spins, and 
later by Barnescl, who was the first to define and use the 
pseudo particle representation of a local impurity level 
we will be working with below. 

Somewhat later these representations were applied to 
the N-fold degenerate Anderson-impurity models in con- 
junction with a rnean-field approximation (MFA) of the 
slave boson ficldgO'Q. As will be discussed below, the 
MFA, inspite of its initial success, suffers from severe 
problems: (1) it leads invariably to the Fermi liquid 
ground state, even for multi-channel models, in contra- 
diction to the exact non-Fermi liquid behavior in those 
cases; (2) even for the single-channel case the MFA fea- 
tures a spurious first order transition separating the low 
temperature Kondo screened state from the high tem- 
perature local moment regime. On the technical side, the 
MFA violates the gauge symmetry mentioned above. The 
latter problem may be corrected by including the leading 
fluctuations around the saddle-point in 1 /N expansionQ, 
whereas the former problems remain. We shall therefore 
avoid using the mean-field approximation which solves 
the problem of the spurious phase transition. 



II. SINGLE- AND MULTI-CHANNEL QUANTUM 
IMPURITY MODELS 

Quantum impurity models such as the Kondo and An- 
derson impurity models were discussed first in the con- 
text of magnetic impurities in metalsu. In these models 
an N-fold spin degenerate local impurity level (here called 
d-level) hybridizes with the conduction band. In the pres- 
ence of a sufficiently strong Coulomb repulsion U between 
any two electrons in the local impurity level and provided 
the level lies sufficiently far below the Fermi energy the 
local level will be singly occupied and a magnetic mo- 
ment forms. The antiferromagnetic exchange interaction 
between the local moment and the local conduction elec- 
tron spin density leads to the Kondo effect. The physics 
of the Kondo effect is not restricted to magnetic ions 
in metals. It depends only on the local level being de- 
generate (any two-level or multi-level system qualifies) 
and the local Coulomb interaction U being sufficiently 
strong. The Kondo effect can therefore occur as well 
for quadrupcde moments, the two-level systems in amor- 
phous. metalsD or for the energy level system of quantum 
dotal3. We will cepsider the prototype Anderson impu- 
rity Hamiltoniantil 



H = Hc + Hd + Hhyb (1) 
with the conduction electron Hamiltonian 

Hc^Y.'>^cj:^^k. ; a = l,...iV (2) 

k,(7 

the local d-electron Hamiltonian 

ndaUda' (3) 

(T a<.<j' 

and the hybridization term 

i^/^yfa = ^^(c^ d.+/i.c) (4) 

ha 

In lowest order perturbation theory in V the d-level at 
energy Ed acquires a width F = 7ry^A/'(0), with A/'(0) 
the conduction electron density of states at the Fermi 
level w = 0. In the Kondo limit of the model, defined by 
Ed < 0, 2Ed + U > 0, F/ I Ed |< 1 and r/{2Ed + U) < 1, 
the low energy sector of the Anderson model may be 
mapped onto the Kondo mode]ll3 featuring a local anti- 
ferromagnetic exchange interaction J of the local spin of 
the single electron in the d-level with the local conduc- 
tion electron spin, where J = V'^[1/ \ Ed\ +l/{2Ed+U)]. 
These models have been studied_extensively by means of 
Wilson'.s xepormalization groupcJ, by the Bethe ansatz 
methodliaij and by means of a phenomenological Fermi 
liquid theoryllEI. Ln this way the following physical pic- 
ture has emergedQ: The model contains a dynamically 
generated low temperature scale, the Kondo tempera- 
ture, which, for [/ — > oo, is expressed in terms of the 
parameters of the Anderson Hamiltonian as 

Tk = DiNT/D)<^^'/^h-^''^/^^, (5) 

with 7V(0) and D = 1/7V(0) the density of states at the 
Fermi energy and the high energy band cutoff, respec- 
tively. A^, M are the degeneracy of the local level and 
the number of conduction electron channels (see below). 

In the intermediate temperature regime, T ^ Tk, res- 
onant spin flip scattering of electrons at the Fermi sur- 
face off the local degenerate level leads to logarithmic 
corrections to the magnetic susceptibility x(^)j the lin- 
ear specific heat coefficient 7(T) and the resistivity p(T) 
and to a breakdown of perturbation theory at T ~ Tk- 
In the single-channel case (M — 1), below Tk a col- 
lective many-body spin singlet state develops in which 
the impurity spin is screened by the conduction electron 
spins as lower and lower energy scales are successively 
approached, leaving the system with a pure potential 
scattering center. The spin singlet formation is sketched 
in Fig. y a) and corresponds to a vanishing entropy at 
T = 0, S{0) ~ 0. It also leads to saturated behavior of 
physical quantities below Tk, like x{T) = const,j{T) = 
c{T)/T = const, and Ap{T) = p{T) - p(0) cx T^, i.e. to 
Fermi liquid behavior. 
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FIG. 1. Sketch of the renormaUzation group for a) the single-channel Kondo model (local moment compensation) and b) 
the two-channel Kondo model (local moment over-compensation). Small arrows denote conduction electron spins 1/2, a heavy 
arrow a localized spin 1/2. The curved arrows indicate successive renormalization steps. 



Multichannel quantum impurity models are character- 
ized by the completely symmetric coupling of M degen- 
erate conduction bands to the impurity. Let us assume 
for definiteness that the ground state of the impurity is a 
magnetic singlet, with M-fold orbital degeneracy, | >, 
where ^ — 1, . . . A/, and the first excited state, obtained 
by adding one electron to the impurity, is magnetic and 
N-fold degenerate, but is an orbital singlet, | T'a >, 
a = l,...,N. An SU(N) x SU(M) symmetric Andetr 
son model may then be defined in the limit C/ — s- oo asB 

H = H,^^ +Ed J2 I {T'a I +Hhyb (6) 

(7 

where 

Hcu^^ezcl ct (7) 

describes the conduction electron continuum in orbital 
channel p, and 

Hhy, ^VY, (c+ ^ I T^i) {T'a \ +h.c) (8) 

k.{T,^ 

and the effect of the Coulomb interaction is to restrict 
the average occupation to Ud = X^o- I r'cr)(r'cr |< 1. 

Like in the single channel case, in the multi-channel 
case, too, the conduction electrons of each channel sepa- 
rately screen the impurity moment by multiple spin scat- 
tering at temperatures below the Kondo scale Tk- How- 
ever, in this case, the local moment is over-compensated, 
since the impurity spin can never form a singlet state 
with both conduction electron channels at the same time 
in this way, as can be seen in Fig. |l| b). As a conse- 
quence of this frustration, there is not a . uiiiq ue ground 
state, leading to a finite residual entropyllZrIla at T = 
of S{Q) = fcslnv^ in the two-channel model. In partic- 
ular, the precondition of FL theory of a 1:1 correspon- 
dence between interacting and non-interacting states is 



violated. As a conseoueace, characteristic singular tem- 
perature dependencetO'tJ of physical quantities persists 

for T ^ Tk down to T = 0: xiT) oc -ln{T/TK), 
c{T)/T « -HT/Tk) and p{T) - p(0) cx -^/tJT^. 



III. PSEUDO-PARTICLE REPRESENTATION 

As discussed above, the dynamics of an electron oc- 
cupying a local level of a quantum impurity will depend 
in an essential way on whether the level is singly oc- 
cupied or multiply occupied, provided the Coulomb in- 
teraction U between two electrons in the same level is 
large. It is therefore useful to divide the Hilbert space 
into sectors labeled by the occupation number of the 
quantum impurity. For each Fock state | a > of the 
impurity one may define a creation operator ajj^, which, 
when operating on a vacuum state | vac >, creates 
the state | a > : | a >= a+ | vac >. For a 
usual impurity level with spin degeneracy N = 2 one 
thus defines two bosons b and a, and two fermions /g-, 
a =tj ii creating the empty and doubly occupied states 
I 6+ I vac >, I 2 0+ I vac >, and the singly 
occupied states | cr >= /+ | vac >. The creation of an 
electron in the empty level, effected by a fermion cre- 
ation operator djr, according to | cr >= d+ | > or 
I 2 >= |i>= —d^ |T>i is described in terms of the 
pseudo-particle operators as 

4 =/+5 + r;„a+/_, (9) 

where the factor rj^ = ±1 for a —] ^ [ accounts for the an- 
tisymmetry of the doubly occupied state. By the intro- 
duction of the pseudo-particles the Fock space has been 
vastly extended. The physical sector of this artificial Fock 
space is defined by the requirement that the impurity 
should be occupied by exactly one pseudo-particle at any 
time (corresponding to a definite state), as expressed by 
the constraint 



fj^fa + b+b 



(10) 



{A)G=tr[pGA]. 



(14) 



In the following we will largely consider the case of in- 
finite Coulomb repulsion [/, implying that double (mul- 
tiple) occupancy of the impurity is not possible. The 
operators a will therefore not appear in this case. 
In the multi-channel SU(N) x SU(M) Anderson model 
discussed in the last section, we will, however, need to 
introduce M boson operators /i — 1, . . . M creating 
the M ground states | F/i >, and N fermion operators 
/+, cr = 1, . . . TV. In terms of these operators the Hamil- 
tonian of the SU(N)xSU(M) Anderson model takes the 
form 



k,cr,tJ. 



(11) 



In order for H to be SU(M) invariant, the slave boson 
multiplet bfi transforms according to the conjugate rep- 
resentation of the SU(M). In addition, the operator con- 
straint 



(12) 



has to be satisfied at all times. One might interpret the 
constraint as a statement of charge quantization, with 
the integer Q the conserved, quantized charge. Simi- 
lar to quantum field theories with conserved charges, the 
charge conservation is intimately related to the existence 
of a local gauge symmetry Indeed, the system defined 
by the Hamiltonian Eq. ( |ll|) is invariant under simulta- 
neous local C/(l) gauge transformations fa fa^'^'^^'^K 
bp, &pe*'^('^', with </i(t) an arbitrary time dependent 
phase. 



A. Exact projection onto the physical Hilbert space 

While the gauge symmetry guarantees the conserva- 
tion of the quantized charge Q, it does not single out 
any particular Q, such as Q = 1. In order to effect the 
projection onto the sector of Fock space with Q =ji 1, one 
may use a procedure first proposed by AbrikosovD: Con- 
sider first the grand-canonical ensemble with respect to 
Q, defined by the statistical operator 



1 



PG 



(13) 



where Zq = tr[exp{— /3(7J-|-A(5)}] is the grand-canonical 
partition function with respect to Q, — A is the associated 
chemical potential, and the trace extends over the com- 
plete Fock space, including summation over Q. The ex- 
pectation value of an observable A in the grand-canonical 
ensemble is given by 



The physical expectation value of A, (A), is to be eval- 
uated in the canonical ensemble where Q = 1. It can 
be calculated from the grand-canonical ensemble by dif- 
ferentiating with respect to the fugacity ( — e~^^ and 
taking A to infinitjCl, 



lim 

A — >oo 



^tr[ie-'5(^+^«)] 
^tr[e-'3(^+^Q)] 



= lim 



{QA)c 

(Q) 



G 



(15) 



Projecting operators acting on the impurity states: — 
We list two important results, which follow straightfor- 
wardly from Eq. (p^): First, the canonical partition func- 
tion in the subspace Q = 1 is 



Zc — lim ti'lQ 

A— ^oo 



-l3iH+\{Q-l))l 



lim (e^^(Q)G(A))ZQ^o , 

A— ^oo 



(16) 



where the subscripts G and C denote the grand-canonical 
and the canonical {Q = 1) expectation value, respec- 
tively. Second, the canonical Q = 1 expectation value of 
any operator A which has a zero expectation value in the 
(5 = subspace, A\Q = 0) — 0, is given by. 



c 



lim 



{A)gW 

{Q)gW 



(17) 



Note that A\Q = 0) — holds true for any physically 
observable operator acting on the impurity. Examples 
are the physical electron operator dj,^ = fabp, or the lo- 
cal spin operator S = X^o-dr' ^fa''^crcr'fcr'- In this case the 
operator Q appearing in the numerator of Eq. ( |l5|) is not 
necessary to project away the Q = sector. In particu- 
lar, the constrained d-electron Green's function is given 
in terms of the grand-canonical one {Gdit-^, T, A)) as 



Gdjuj^T, A) 



(18) 



In the enlarged Hilbert space {Q — 0, 1, 2, ...) Gd{uj, T, A) 
may be expressed in terms of the grand-canonical 
pseudo-fermion and slave boson Green's functions using 
Wick's theorem. These auxiliary particle Green's func- 
tions, which constitute the basic building blocks of the 
theory, are defined in imaginary time representation as 

GfaiTl - T2) = -(f {/.(Ti)/t(r2)}>G (19) 

Gbpiri - r2) = -{f{b4n)blir,)})G , (20) 

where T is the time ordering operator. The Fourier trans- 
forms of Q f i, may be expressed in terms of the exact self- 
energies T,f ii as 

GfA^^n) = [[GltiiL^n]]-^ - E/,b(ic^„)}"' (21) 



4 



where 



Ed-\y 



(22) 
(23) 



Since, as a consequence of the projection procedure A 
oo , the energy eigenvalues of H + XQ scale to infinity as 
AQ, it is useful to shift the zero of the auxiliary particle 
frequency scale by A (in the Q = 1 sector) and to define 
the "projected" Green's functions as 



Gf^b{^) = lim Qffiii^ 

A — ^oo 



A) 



(24) 



Note that this does not affect the energy scale of physi- 
cal quantities (like the local d electron Green's function), 
which is the difference between the the pseudo-fermion 
and the slave-boson energy. 

Canonical expectation values of conduction electron op- 
erators: — The canonical (i.e. projected onto the Q — 1 
subspace), local conduction electron Green function is 
given as 



Gc^^ai^^^n) = \[G°^„{iUJn)] ^-ScMa(«^n)| (25) 



with 



k k 

(26) 



where /Zc is the chemical potential of the conduction 
electrons. The canonical, local c-electron self-energy, 
Sc^cr(iw„), cannot be obtained from the grand-canonical 
one by simply taking the limit A — > oo, since the c- 
electron density has a non-vanishing expectation value 
in the Q = subspace. However, it follows immediately 
from the Anderson Hamiltonian, Eqs. (§)-(§), that the 
exact, canonical conduction electron t-matrix t^^{iuj), 
defined by Gc — G"[l + tG"], is proportional to the 
full, projected c?-electron Green's function, t^„{iuj) = 
\V\^Gd^rj- Thus, we have as an exact relation, 



Gc^,a{iuJn) = G° (iw„) 1 + |l/|2Gd^.(iw„)G° (iw„) 



(27) 



and by comparison with Eq. ( p5[ ) we obtain the local con- 
duction electron self-energy respecting the constrained 
dynamics in the impurity orbital. 



\V?Gd^„{iUr.) 



l + \V\^G%,{^uJ^)Gd^.a{^iOr,) 



(28) 



Using phenomenological Fermi liquid theoryll3 and also 
by means of pertuijbation theory to infinite order in the 
on-site repulsion UYn it has been shown for the Fermi liq- 
uid case M = 1 of the symmetric Anderson model {2Ed = 
—U) that the exact d-electron propagator Gda{'^) and 
the d-electron self-energy ^^^(a;) = to — Gdcr{t^)~^ obey 
the following local Fermi liquid relations in the limit 



Luttinger theorem: 



duf{uf-^f^Gd.{oj)^G 



Friedel-Langreth: — ImG(jcr(w) = 4 sin^ \ — c \( + ('!^\ 

TT 



Im^dai^j) 



V \ N J 
T 



sin\-Knd/N) \iiin\-Knd/ N) 



\TkJ 
T 



2 



T, 



K 



(29) 
(30) 
(31) 



where c is a constant of 0(1). Combining Eqs. (28), 
(|^) it follows that (for M = 1, away from particle 
hole symmetry) Sco- exhibits (in an exact theory) local 
Fermi liquid behavior as well, lmTjca{^^ — iO, T = 0) = 
a + 6(a;/Tfj-)^ for w — > 0. Note that this quantity is 
different from the grand-canonical conduction electron 



self-energy and has a finite imaginary part at the Fermi 
level. 

The momentum dependent conduction electron 
Green's function in the presence of a single impurity 
is given in terms of the canonical d-electron propagator 
as 



Gc^,a{k,k' ] icj„)G° (fc,iw„) % -f |ypGd(it^„)G° (fc', itj„) 



(32) 



The latter expression is the starting point for treating a random system of many Anderson impuritiesi 



i. 
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B. Analytical properties and infrared behavior 

The Green's functions Gf^b,c have the following spec- 
tral representations 

Gf,b,Mu;„) = r du'^^^^^ (33) 

with the normalization of the spectral functions ^/,6,c 

/oc 
dujAf^kA^) = 1. (34) 
-oo 

Taking the limit \ ^ oo has important consequences on 
the analytical structure of the auxiliary particle Green's 
functions: 



(1) — It follows directly from the definitions Eqs. (|24|), 
(^), using Eqs. (|l3|), (|lj), that the traces appearing in 
the canonical functions G/A-.are taken purely over the 
(5 = sector of Fock spaca£3. Thus, the backward-in- 
time (ti < T2) or hole-like contribution to the auxiliary 
propagators in Eq. (^^ vanishes after projection, and 
we have 

G/.(ti - T2) = -e(ri - T2) lim (/,(Ti)/t(T2))G (35) 

A—* 00 

Gf,^(ri - T2) = -e(ri - T2) lim (6p(ri)6t (r2))G • (36) 

Consequently, their spectral functions Af^i, have the 
Lehmann representation 

(37) 
(38) 



I 

= ^ e-P^''-\{lMfl\Q.rn)\^5{u-{El-El)) 

m,n>0 

= e-'^^-|(l,n|6t|0,m)p5(^-(£;i-O), 

m,n>0 



where E'^ are the energy eigenvalues (iJg < £','3 is the 
ground state energy) and \Q,n) the many-body eigen- 
states of H in the sector Q of Fock space. At zero temper- 
ature, Af reduces to Af^{ijj) = J2n>o I (1, '^l/J|0, 0) p 
S{uj — {E^ — Eq)) and similar for Ai,. It is seen that the 
Af^b have threshold behavior at lu = Eq = Eq — Eq, with 
Af^h{uj) = for CO < Eq, T = 0. The vanishing imag- 
inary part at frequencies w < may be shown to be a 
general property of all quantities involving slave particle 
operators, e.g. also of auxiliary particle self-energies and 
vertex functions. 

(2) — As will be seen in section V B (Eq. (|5^)), phys- 

Eq is the ground state energy in the Q = I sector. The 
expressions (|37|), ( |3^ ) and (pl|), ( ^2| ) immediately imply 
a relation between Af^b and Aj^^, 

Aj.iio) = e~P-Af,b{io) . (43) 

(3) — The prop erty of only forward-in-time propa- 
gation (Eqs. (p5|), (pq)) means that the auxiliary parti- 
cle propagators Gf^t are formally identical to the core 
hole prooagators of the well-known X-ray threshold 
problemE3-c3. Thus, the knowledge of the infrared be- 
havior of the latter may be directly applied to the for- 
mer. In particular, the spectral functions are found (see 
below) to diverge at the threshold Eq in a power law 



ical expectation values not only involve the particle-like 
auxiliary propagators Eqs. (p5|), ( ^6| ) but also hole-like 
contributions. It is, therefore, useful to define the "anti- 
fermion" and "anti-boson" propagators (in imaginary 
time representation) 

GJ,{t, - r2) - -e(r2 - n) hm (/] (t2)/.(ti))g (39) 
G^-^n - T2) = -e(r2 - n) hm (4(r2)6^(ri))G , (40) 

whose spectral functions have the Lehmann representa- 
tions 

(K-ElJ) (41) 
iK-E'J). (42) 

I 

fashion (infrared singularity) 

Af^i,{uj)^\iu-EQ\-"f-'0{Lu-Eo) (44) 

due to a diverging number of particle-hole excitation pro- 
cesses in the conduction electron sea as oj Eg- For the 
single channel case (M = 1), i.e. the usual Kondo or 
mixed valence problem, the exponents a/ and at can 
be found analytically from the following chain of argu- 
ments: Anticipating that in this case the impurity spin 
is completely screened by the conduction electrons at 
temperature T = 0, leaving a pure-potential scattering 
center, the ground state |1,0) is a Slater determinant of 
one-particle scattering states, characterized by scattering 



I 

m,n>0 
m,n>0 
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phase shifts ry^ in the s-wave channel (assuming for sim- 
phcity a momentum independent hybridization matrix 
element V). To calculate the fermion spectral function 
Afcr{uj) at T = from Eq. (|3^), one needs to evaluate 
(1, n \ fl \ 0,0), which is just the overlap of two slater de- 
terminants, an eigenstate of the fully interacting Kondo 
system, on the one hand, and the ground state 

of the conduction electron system in the absence of the 
impurity combined with the decoupled impurity level oc- 
cupied by an electron with spiiincr, /||0,0), on the other 
hand. As shown by AndersorO, the overlap of the two 
ground state slater determinants, (1,0 | | 0,0), tends 
to zero in the thermodynamic limit (orthogonality catas- 
trophe). Analogous relations hold for the boson spectral 
function Ab{uj). As a result, the exponential relaxation 
into the interacting ground state for long times is inhib- 
ited, leading to the infrared power law divergence of the 
spectral functions, Eq. (^). 

The X-ray threshold exponents can be expressed in 
terms of the scattedng phase shifts at the Fermi level by 
the exact relationEEl 



1 X^Cnf.ba' 



(45) 



using the Bethe|-ansatz solution and boundary confor- 
mal field theoryl^. It is interesting to note that (i) 
the exponents depend on the level occupancy (in the 
Kondo limit Ud ^ 1, aj = 1/N and — I — ^/N, 
whereas in the opposite, empty orbital, limit Ud — > 0, 
a/ — > and a?, — )■ 1) (ii) the sum of the exponents 



1 



)nd(l-nd) 



> 1. 



af + ab — ^ I ^ 

We stress that the above derivation of the infrared ex- 
ponents a-f^b holds true only if the impurity complex acts 
as a pure potential scattering center at T = 0. This 
is equivalent to the statement that the conduction elec- 
trons behave locally, i.e. at the impurity site, like a Fermi 
liquid. Conversely, in the multi-channel (non-FL) case, 
N > 2, M > N, the exponents-have been found from a 
conformal field theory solutiorO of the problem in the 
Kondo limit to be 



af = M/{M + N) 
ab = N/{M + N) 



(50) 



which differ from the FL values. Thus, one may infer 
from the values of a/^f, as a function of Ud, whether or 
not the system is in a local Fermi liquid state. 



Here the 77/0-' (?7f)o-') are the scattering phase shifts of 
the single-particle wave functions in channel a' of the 
fully interacting ground state |1,0), relative to the wave 
functions of the free state /]|0,0) (61'|0,0)). Via the 
Friedel sum rule, the scattering phase shifts are, in turn, 
related to the change Ariay' of the average number of 
conduction electrons per scattering channel a' due to the 
presence of the impurity: r]f,ba' = T^^Tica'- Obviously, 
Ariccr' is equal and opposite in sign to the difference of 
the average impurity occupation numbers of the states 
|1,0) and /J|0,0) (6^|0,0) ). Thus, in the pseudofermion 
propagator Gfa we have the phase shifts. 



and in the slave boson propagator Gb, 

nd 

Vb = , 



(46) 



(47) 



where Ud denotes the total occupation number of the im- 
purity level in the interacting ground state. (The term 
6cra' in ?7/ cr' appcars because /j^|0,0) has impurity occu- 
pation number 1 in the channel a.) For example, in the 
Kondo limit ^ 1 and for a spin 1/2 impurity (A^ — 2) 
this leads to reamance scattering, rjf ba' = 7i'/2. As a 
result, one find£3 for the threshold exponents 



ab 



2nd 



1 



N 
N 



(48) 
(49) 



These results have been found independently fii©i«-|Wil- 
son's numerical renormalization group approachc§E3 and 



IV. 



MEAN FIELD APPROACH AND 1/N 
EXPANSION AT (7 ^ 00 



d hybridiza- 



tion of the Anderson model (|il|) is much smaller than the 

where AA(0) = l/D 



For physical situations of interest, the s - 

(p]) is mT 

conduction band width, M{0)V <C 1, 
is the local conduction electron density of states at the 
Fermi level. This suggests a perturbation expansion in 
N{0)V . A straightforward expansion in terms of bare 
Green's functions is not adequate, as it would not allow 
to capture the physics of the Kondo screened state, or else 
the infrared divergencies of the auxiliary particle spectral 
functions discussed in the last section. In the framework 
of the slave boson representation, two types of nonper- 
turbative approaches have been developed. The first one 
is mean field theory for both the slave boson amplitude 
(6) and the constraint ((Q) = 1 rather than Q = 1). The 
second one is resummation of the perturbation theory to 
infinite order. 



A. Slave boson mean field theory 

Slave boson mean field theory is based on the assump- 
tion that the slave bosons condense at low temperatures 
such that {bp) ^IL Replacing the operator 6^ in H + \Q 
by (see Ref.El), where A is a Lagrange multiplier to 
be adjusted such that (Q) = 1, one arrives at a reso- 
nance level model for the pseudofermions. The position 
of the resonance, Ed -I- A, is found to be given by the 
Kondo temperature Tk, and is thus close to the Fermi 
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energy. The resonance generates the low energy scale 
Tk, and leads to local Fermi liquid behavior. While this 
is qualitatively correct in the single-channel case, it is 
in blatant disagreement with the exactly known behav- 
ior in the multi-channel case. The mean field theory can 
be shown to be exact for M = 1 in the limit N ^ oo 
for a model in which the constraint is softened to be 
Q = N/2. However, for finite N the breaking of the local 
gauge symmetry, which would be implied by the conden- 
sation o£J;he slave boson field, is forbidden by Elitzur's 
theorems. It is known that for finite N the fluctuations 
in the phase of the complex expectation value (bp) are 
divQwfcmt and lead to the suppression of (bp) to zero (see 
alsoE^-Ea). This is true in the cartesian gauge, whereas in 
the radial gauge the phase fluctuations may be shown to 
cancel at least in lowest order. It has not been possible 
to connect the mean field solution, an apparently rea- 
sonable description at low temperatures and for M = 1, 
to the high temperature behavior (T 3> TV), dominated 
by logarithmic temperature dependence, in a continuous 
wa5iO. Therefore, it seems that the slave boson mean 
field solution does not offer a good starting point even 
for only a qualitatively correct description of quantum 
impurity models. 

B. 1/N expansion vs. self— consistent formulation 

The critical judgement of mean field theory is corrob- 
orated by the results of a straightforward 1 /A^-expansion 
in the single channel case, keeping the exact constraint, 
and |Bot allowing for a finite bose field expectation 
valueEa. Within this scheme the exact behavior of the 
thermodynamic quantities (known from the Bethe ansatz 
solution) at low temperatures as well as high tempera- 
tures is recovered to the considered order in Also, 
the exact auxiliary particle exponents a/^b are repro- 
duced m, order 1/iV, using a plausible exponentiation 
schemeEZl. 

In addition, dynamical quantities like the d-electron 
spectral function and transport coefficients can be cal- 
culated exactly to a desired order in within this 
approach. However, as clear-cut and economical this 
method may be, it does have serious limitations. For 
once, the 1/7V expansion is^ot uniformly convergent as 
a function of temperatureE3. Rather, the expansions 
at low temperature and at high temperature have to 
be done around two different saddle-points (the limits 
N —f oo). It is not known how to match these expan- 
sions in the crossover region T ~ Tk in a systematic 
way. Secondly, the experimentally most relevant case of 
A'' = 2 or somewhat larger is not accessible m 1/N ex- 
pansion. Thirdly, non-Fermi liquid behavior, being nec- 
essarily non-perturbative in cannot be dealt with 
in a controlled way on the basis of a 1/A'^-expansion. 
To access these latter two regimes, a new approach non- 
perturbative in 1/N is necessary. 



V. CONSERVING APPROXIMATIONS: GAUGE 
INVARIANT SELF-CONSISTENT 
PERTURBATION THEORY IN THE 
HYBRIDIZATION 



We conjecture that this new approach is gauge in- 
variant many-body theory of pseudofermions and slave 
bosons. As long as gauge symmetry violating objects 
such as Bose field expectation values or fermion pair cor- 
relation functions do not appear in the theory, gauge 
invariance of physical quantities can be guaranteed in 
suitably chosen approximations by the proper match of 
pseudofermion and slave boson properties, without intro- 
ducing an additional gauge-field. This requires the use of 
conserving approximationsLJ, derived from a Luttinger- 
Ward functional <I>. 



A. Generating functional 

$ consists of all vacuum skeleton diagrams built out of 
fully renormalized Green's functions Gb.f,c and the bare 
vertex V. The self-energies 'Sbj,c are obtained by taking 
the functional derivative of $ with respect to the corre- 
sponding Green's function (cutting the Green's function 
line in each diagram in all possible ways), 

Sh,/,c = S^/SGbj,c- (51) 

Irreducible vertex functions, figuring as integral kernels 
in two-particle Bethe-Salpeter equations, are generated 
by second order derivatives of 

The choice of diagrams for $ defines a given approx- 
imation. It should be dictated by the dominant physi- 
cal processes and by expansion in a small parameter, if 
available. Even in the presence of a small parameter (in 
this case AfaV), straightforward low order renormalized 
perturbation theory may not give an even qualitatively 
correct result, if singular vertex functions appear. It is 
therefore necessary to check whether vertex functions be- 
come singular at the level where approximate single par- 
ticle Green's functions are used to calculate the integral 
kernels of the respective vertex functions. Should this be 
the case, the vertex functions have to be included into the 
approximation in a gauge invariant way. As wc shall see, 
in the present case it is required to take all two-particle 
vertex functions into consideration. Vertex functions in- 
volving three or more interacting particles will be omitted 
in the hope that the corresponding phase space is small 
so that they will contribute less even if they are singular. 
This is the first fundamental approximation. The second 
one is that we will approximate the irreducible kernels in 
the Bethe-Salpeter equations for the vertex functions by 
the lowest order diagram in NqV. 
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FIG. 2. Diagrammatic representation of the generating functional <1> of the NCA. Also shown are the pseudoparticle 
self-energies and the local electron Green's function derived from $, Eqs. (19)-(21). Throughout this article, dashed, wavy 
and solid lines represent fermion, boson, and conduction electron lines, respectively. In the diagram for E /o- the spin labels are 
shown explicitly to demonstrate that there are no coherent spin fluctuations taken into account. 



B. Non-crossing approximation (NCA) 

As noted before, in the present context, we may take 
the hybridization y to be a small quantity (dimensionless 
parameter NqV). This suggests to start with the lowest 
I 



where = i ImG'2^^/A/'(0) is the (free) conduction 

electron density of states per spin and channel, normal- 
ized to the density of states at the Fermi level A/'(0), 
and f{e) = l/(exp(/3£) + 1) denotes the Fermi distribu- 
tion function. Together with the expressions (pi]), ( |2^ ) 
for the Green's functions, Eqs. (p2|)-(p4|) form a set of 
self-consistent equations for Sbj.c, comprised of all di- 
agrams without any crossing propagator lines and are, 
thus, knovHipa^ the "non-crossing approximation" , in 
short NCAm. 

At zero temperature and for low frequencies, Eqs. js^ ) 
and (^3|) may be converted iiitp a set of linear differen- 
tial equations for G/ and G(,c3, which allow to find the 
infrared exponents as a/ = m+n ' '^^ ~ ai+n ' indepen- 
dent of Ud- For the single channel case these exponents 
do not agree with the exact exponents derived in section 
III B. This indicates that the NCA is not capable of re- 
covering the local Fermi liquid behavior for M = 1. A 
numerical evaluation of the d-electron Green's function, 
which is given by the local self-energy Ec divided by 
and hence is given by the boson-fermion bubble within 
NCA (Fig. ^), shows indeed a spurious singularity at the 



order (in V) diagram of which is second order (see 
Fig. H) . The self-energies generated from this obey after 
analytic continuation to real frequencies {iuj —^uj — iO) 
and projection the following equations of self-consistent 
second order perturbation theory 



(52) 
(53) 
(54) 



Fermi energyH^. The NCA performs somewhat better in 
the multi-channel case, where the exponents a/ and at 
yield the correct non-Fermi liquid exponents of phyai^ 
cal quantities as known froia the Bethe ansatz solutionEj 
and conformal field theorylij. However, the specific heat 
and the residual entropy are not given correctly in NCA. 
Also, the limiting low temperature scaling laws for the 
thermodynamic quantities are attained only at temper- 
atures substantially below Tk, in disagreement with the 
exact Bethe ansatz solution. 



C. Evaluation of the self—consistency equations at 
low temperatures 

In order to enter the asymptotic power law regime of 
the auxiliary spectral functions, the self-consistent equa- 
tions must be evaluated for temperatures several orders 
of magnitude below Tk, the low temperature scale of the 
model. The equations are solved numerically by itera- 
tion. In the following we describe the two main proce- 
dures to make the diagrammatic auxiliary particle tech- 



Si;^^^)(^-z0)=r^ J ^/(£)A^^^(£)G/.(a; + £-zO) 

G^a^^'' - *0) = / dee-P'[GfA^ + e - iO)A,,(e) - Af„{e)Gbf.{e - + iO)] 
^ Jde [GMu; + e - iQ)A--{e) - Aj^{e)Gtf.{s - u; + lO)] , 
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nique suitable for the lowest temperatures. 

The grand-canonical expectation value of the auxiliary 
particle number appearing in Eq. (]l^) is given in terms 
of the grand-canonical (unprojected) auxiliary particle 
spectral functions Af^bi^, A) by, 

(Q)g(A) = (55) 



duj 



Jiuo) J2 Af. [uo, A) + h{Lo) J2 Am A)' 



where /(w), &(w) are the Fermi and Bose distribution 
functions, respectively. Substituting this into the expres- 
sion ( [l6| ) for the canonical partition function we obtain 
after carrying out the transformation — > + A, and 
taking the limit X oo 



-/3F.™j,(T) 



Zq=o 



lim el^^{Q)G{X) 



(56) 



By definition Fi^p — — ^ ln(Zc/^Q=o) is the impurity 
contribution to the Free energy. 

The numerical evaluation of expectation values like 
(g)G(A ^ oo) (Eq. @)) or GdM,(a;, A oo) (Eq. H)) 
is non-trivial, (1) because at T = the auxiliary spec- 
tral functions Af^b{uj,T) are divergent at the threshold 
frequency Eq, where the exact position of Eq is a priori 
not known, and (2) because the Boltzmann factors e~^'^ 
diverge exponentially for to < 0. Therefore, we apply the 
following transformations: 

(1) Before performing the projection u uj + X, 
X oo we re-define the frequency scale of all auxil- 
iary particle functions Af^b according to w uj + Xq, 
where Aq is a finite parameter. In each iteration Aq is 
then determined such that 



<k.e-P-[Y,Af^{Lo)+J2Ab,ico) 



1 



(57) 



where Af^^i^) = hmA^oo Af,i,{uj -I- Aq, A) is now an aux- 
iliary spectral function with the new reference energy. 



It is seen by comparison with Eq. ( pq ) that Ao(T) = 
F.mpiT) = FQ=i(r) - Fq^o{T), i.e. Aq is the chem- 
ical potential for the auxiliary particle number Q, or 
equivalently the impurity contribution to the Free energy. 
The difference of the Free energies becomes equal to the 
threshold energy Eq = E^^^ - E'gfp at T = 0. More 
importantly, however, the above way of determining a 
"threshold" is less ad hoc than, for example, defining it 
by a maximum in some function appearing in the NCA 
equations. It is also seen from Eq. ( ^ ) that this pro- 
cedure defines the frequency scale of the auxiliary par- 
ticles such that the T = threshold divergence of the 
spectral functions is at the fixed frequency w = 0. This 
substantially increases the precision as well as the speed 
of numerical evaluations. Eq. (18) for the projected d- 
electron Green's function becomes 



Gdico) = lim e>'^GdiLu,T,X). 

X—>-oo 



(58) 



(2) The divergence of the Boltzmann factors implies 
that the self-consistent solutions for A/j,(cj) vanish ex- 



ponentially 



for negative frequencies, confirming 



their threshold behavior. It is convenient, not to formu- 
late the self-consistent equations in terms of Af^i, like in 
earlier evaluationsa, but to define new functions Af,i,{oj) 
and ImSy.b(a;) such that 



ImE/,f,(Lj) = f{-uj) ImE/,6(cj). 



(59) 
(60) 



After fixing the chemical potential Aq and perform- 
ing the projection onto the physical subspace, the 
canonical partition function (Eq. (^6|)) behaves as 
limA^oo e^'-^-^"'^ Zc{T) = 1, and from Eq. (|43|) we have 
A^ ^{lo) = f{ijj)Af,b{ijj)- In this way all exponential di- 
vergencies are absorbed by one single function for each 
particle species. As an example, the NCA equations in 
terms of these functions are free of divergencies of the 
statistical factors and read 



f[-e){l-f{u~s)) , ~ 

— AcaA^)Abf,{UJ - £) 



/(£)(l-/(c. + e)) ^0 



ImS,p(a;-»0)=rX] 
(Q)(Ao,A^^)= j dujJ{uj)[Y,Af,{io) + Y.^bf.{' 



A%A^)Af^{u + e) 



= 1 



ImGrf,(c^ -iO)^ / de[f{e + u)f(~e) + /(-e - c^)/(e)]i/,(e + u)Ab{e). 



(61) 
(62) 
(63) 
(64) 
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a |i a |i' co'-^2 c b c c c 




FIG. 3. Diagrammatic representation of the Bethe-Salpeter equation for (1) the conduction electron-pseudofermion 
T- matrix r*'^^^ Eq. (||), and (2) the conduction electron-slave boson T-matrix T*'^''', Eq. T'"''^ is obtained from T'^^^ 

by interchanging f ^ b and c <-> . 



The real parts of the self-energies S/ , S;, are determined 
from ImE^, ImSb through a Kramers-Kroenig relation, 
and the auxiliary functions Afc^iuj) = -^ImS f^{LLj — 

Abpiuj) = ilmEf,^(cj - iO)/[(cj + Aq - iO - ReT,bfi{uj - 

iO)) + lmSi,fi{uj — «0)^] , thus closing the above set of 
equations. 

The method described above allows to solve the NCA 
equations effectively for temperatures down to typically 
T = 10~^Tr-. It may be shown that the procedures 
described above can also be applied to self-consistently 
compute vertex corrections beyond the NCA (see section 
VI), thus avoiding any divergent statistical factors in the 
selfconsistency equations. 

VI. CONSERVING T-MATRIX 
APPROXIMATION (CTMA) AT f/ ^ oo 

A. Dominant contributions at low energy 

In order to eliminate the shortcomings of the NCA 
mentioned above, the guiding principle should be to find 
contributions to the vertex functions which renormalize 
the auxiliary particle threshold exponents to their correct 
values, since this is a necessary condition for the descrip- 
tion of FL and non-FL behavior, as discussed in section 
III B. Furthermore, it is instructive to realize that in 



NCA any coherent spin flip and charge transfer processes 
are neglected, as can be seen explicitly from Eqs. js^), 
( p3| ) or from Fig. |2| These processes are known to be re- 
sponsible for the quantum coherent collective behavior of 
the Anderson impurity complex below Tk ■ The existence 
of collective excitations in general is reflected in a singular 
behavior of the corresponding two-particle vertex func- 
tions. In view of the tendency of Kondo systems to form a 
collective spin singlet state, we focus our attention on the 
two-particle vertex functions, in particular, in the spin 
singlet channel of the pseudofermion-conduction electron 
vertex function and in the slave boson-conduction elec- 
tron vertex function. It may be shown by power count- 
ing arguments (compare appendix A) that there are no 
corrections to the NCjA exponents in any finite order of 
perturbation theorycS. Thus, we are led to search for sin- 
gularities in the aforementioned vertex functions arising 
from an infinite resummation of terms. 

From the preceding discussion it is natural to perform 
a partial resummation of those terms which, at each or- 
der in the hybridization V, contain the maximum number 
of spin flip or charge fluctuation processes, respectively. 
This amounts to calculating the conduction electron- 
pseudofermion vertex function in the "ladder" approx- 
imation defined in Fig. |^, where the irreducible vertex is 
given by V^Gb- In analogy to similar resummations for 
an interacting one-component Fermi system, we call the 
total c-f vertex function T-matrix T^'^f\ The Bethe- 
Salpeter equation for T^^f^ reads (Fig. | (1)), 



I 

-V^TY,Gb-,{^u^n + ^^'i~^^n)Gf.{^u:':^) Gl^^iz^l^ - iu:'l) r^^^,^:, (zc.;', zc.;, zO„). (65) 
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A similar integral equation holds for the charge fluctuation T-matrix T*^*^^^ (Fig. || (2)), 

(jj'j 



In the above Bethe-Salpeter equations cr, r, cr', r' repre- 
sent spin and /j,, v, v' channel indices. We note that 
these are the only two-particle vertex functions after pro- 
jection. The principal approximation adopted here is the 
form of the irreducible kernel, which we approximate by 
the lowest order diagram. 

Inserting NCA Green's functions for the intermediate 
state propagators of Eq. ( |6^ ) and solving it numerically, 
we find at low temperatures and in the Kondo regime 
{ud ~ 0.7) a pole of T^'^^^ in the singlet channel (see ap- 
pendix A) as a function of the center-of-mass (COM) 
frequency fi, at a frequency which scales with the Kondo 
temperature, 17 = fie/ — —Tk- This is shown in Fig. ^. 
The threshold behavior of the imaginary part of T^^f^ as 
a function of with vanishing spectral weight at nega- 
tive frequencies and temperature T = is clearly seen. In 
addition, a very sharp structure appears, whose broaden- 
ing is found to vanish as the temperature tends to zero, 
indicative of a pole in T^'^^^ at the real frequency r^c/, 
i.e. the tendency to form a collective singlet state be- 
tween the conduction electrons and the localized spin. 

4 I ■ ■ y — r-- ' ■ ' 

>* 

° 3 - 




-10 -5 5 10 

^ID (xlO"^) 

FIG. 4. Imaginary part of the conduction elec- 

tron-pseudofermion T-matrix t''^^' as a function of the COM 
frequency Q. for the single-channel case M = 1, N = 2, eval- 
uated by inserting NCA solutions for the intermediate state 
propagators {Ed = -0.67D, P = 0.15-D, T = 4-10"^Tk). The 
contribution from the pole positioned at a negative frequency 
57 = 57c/ — —Tk (compare text) is clearly seen. 

Similarly, the corresponding T-matrix T'^'^^'^ in the con- 



I 

duction electron-slave boson channel, evaluated within 
the analogous approximation, develops a pole at nega- 
tive values of f2 in the empty orbital regime (n^ ^ 0.3). 
In the mixed valence regime (n^ ~ 0.5) the poles in both 
T'^^^'^ and T^^^) coexist. The appearance of poles in the 
two-particle vertex functions r^^"^^^ and T^'^''\ which sig- 
nals the formation of collective states, may be expected 
to influence the behavior of the system in a major way. 



B. Self— consistent formulation: CTMA 



On the level of approximation considered so far, the 
description is not yet consistent: In the limit of zero 
temperature the spectral weig ht of T('=/) and T'-"''^ at 
negative frequencies f2 should be strictly zero (threshold 
property). Nonvanishing spectral weight at 17 < like a 
pole contribution for negative fl in T^'^^^ or T^'^^) would 
lead to a diverging contribution to the self-energy, which 
is unphysical. However, recall that a minimum require- 
ment on the approximation used is the preservation of 
gauge symmetry. This requirement is not met when the 
integral kernel of the T-matrix equation is approximated 
by the NCA result. Rather, the approximation should be 
generated from a Luttinger-Ward functional. The corre- 
sponding generating functional is shown in Fig. ^ It is 
defined as the infinite series of all vacuum skeleton dia- 
grams which consist of a single ring of auxiliary particle 
propagators, where each conduction electron line spans at 
most two hybridization vertices. As shown in appendix A 
by means of a cancellation theorem, the CTMA includes, 
at any given loop order, all infrared singular contribu- 
tions to leading and subleading order in the frequency uj. 
The first diagram of the infinite series of CTMA terms 
corresponds to NCA (Fig. ||). The diagram containing 
two boson lines is excluded, since it is not a skeleton. Al- 
though the spirit of the present theory is different from 
a large N expansion, it should be noted that the sum 
of the $ diagrams containing up to four boson lipes in- 
cludes all terms of a 1/N expansion up to 0(1/A^^)llj. By 
functional differentiation with respect to the conduction 
electron Green's function and the pseudofermion or the 
slave boson propagator, respectively, the shown $ func- 
tional generates the ladder approximations T^'^^\ T^^^^ 
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NCA 




spin 






charge 



fluctuations 



FIG. 5. Diagrammatic representation of the Lut- 

tinger-Ward functional generating the conserving T-matrix 
approximation (CTMA). The terms with the conduction elec- 
tron lines running clockwise (labeled "spin fluctuations" ) gen- 
erate the T-matrix T^'^-^\ while the terms with the conduc- 
tion electron lines running counter-clockwise (labeled "charge 
fluctuations") generate the T-matrix t''^''). 



for the total conduction electron-pseudofermion ver- 
tex function and for the total conduction electron-slave 
boson vertex function (Fig. ^). The auxiliary parti- 
cle self-energies are obtained in the conserving scheme 
as the functional derivatives of $ with respect to G/ 
or Gb, respectively (Eq. (|l|)). This defines a set of 
self-consistency equations, vifhich we term conserving T- 
matrix approximation (CTMA), where the self-energies 
are given as nonlinear and nonlocal (in time) function- 
als of the Green's functions, while the Green's functions 
are in turn expressed in terms of the self-energies. The 
solution of these equations requires that the T-matrices 
have vanishing spectral weight at negative COM frequen- 
cies fl. Indeed, the numerical evaluation shows that the 
poles of T^'^/) and T^'^^) are shifted to f7 = by self- 
consistency, where they merge with the continuous spec- 
tral weight present for > 0, thus renormalizing the 
threshold exponents of the auxiliary spectral functions. 

In the following we give explicitly the self-consistent 
equations which determine the auxiliary particle self- 
energies within CTMA. For that purpose, it is useful 
to define conduction electron-fermion and conduction 
electron-boson vertex functions T*^'^-^'^*\ (±) with- 
out (-I-) or with (— ) an alternating sign between terms 
with even and odd number of rungs (compare Fig. ^) . In 
the Matsubara representation, these vertex functions, to 
be labelled "even" (-I-) and "odd" (— ) below, are given 
by the following Bethe-Salpeter equations: 



J 



T 



(c/)(±)m 



(67) 



and 

,(cb)(±)a 



(69) 



U'•'''l^^i^u;n,^u;'^,ln^) = -V'^Y^GfA^^n + lu;'; ~ ^n^)Gb4^^':,) G°,,(zc.: - zf}„) G ^lu;'^ + lu;'; - iQ,,) . (70) 



Note that, in addition to the alternating sign, these 
vertex functions differ from the T-matrices defined in 
Eqs. (|65|), ( |66| ) in that they contain only terms with 
two or more rungs, since the inhomogeneous parts [Z'-'^-^' 
and [/('^^J represent terms with two bosonic or fermionic 
rungs, respectively. The terms with a single rung corre- 



I 

spond to the NCA diagrams and are evaluated separately 
(see below). 

The spin degrees of freedom of r('=-'') are uniquely 
determined by the spin indices a, r of the ingoing con- 
duction electron and pseudofermion lines (Fig. 0). It is 
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FIG. 6. Diagrammatic representation of the CTMA expressions for pseudoparticle self-energies E/ and E;,. The first term 
drawn on the righthand side of E/ and E;,, respectively, is the NCA diagram. The diagrammatic parts t['^P , t['^2^ are explained 
in the text, Eqs. (|8|)-(||). 



instructive to note that in the spin 5=1/2 case {N ~ 2) boson vertex function in terms of the channel degrees 
the singlet and triplet vertex functions (which correspond of freedom v. The total CTMA pseudoparticle self- 



to the two-particle Green's functions in the singlet chan- 
nel, 0^ ~ E.(n(c./-.-c_./.)(4/l^-c!.,/t)}), and 
in the triplet channel with magnetic quantum number consist of three terms each, 

TO = o,±i, (/>^^o ~ E.(n(c./-. + c_,/,)(4/L + 

C-<t/])}): ~ (r{c±i/±icl^i/^i}), respectively) 

may be identified in the following way, 



energies, as derived by functional differentiation from the 
generating functional Fig. ^, are shown in Fig. ^ and 
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(c/)t _ 
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(71) 
(72) 
(73) 



Sfc;,(icj„) = E^^'^'^^(iw„) + E[^^^(iw„) + E[^^^(iw„) 



(74) 

)• 

(75) 



The first term of S/ and Sf, represents the NCA self- 
energies, Eqs. (|5^), (|5^). The second and third terms 
arise from the spin and the charge fluctuations, respec- 
tively, and are given for pseudofermions by 



Analogous relations hold for the conduction electron- 
I 



and for slave bosons by 

(iw„) = ivi^Gjl(ic^„-iaOri("''Hjc^„,iw»,if^„) , (79) 



where the vertex functions appearing in these expressions are defined as 

N -1 
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(-) _ ]\,-[ ijicb) 



(82) 

(83) 



These combinations of the even and odd vertex functions 
ensure the proper spin and channel summations in the 
self-energies. For the sake of clarity, the spin and channel 
indices as well as the frequency variables are not shown 
explicitly. In Eqs. (|80|), ( p^ the terms with two rungs, 
N U'^'^f'^ , M , have been subtracted, since they would 
generate non-skeleton self-energy diagrams. Likewise, in 
Eqs. (^l|), ( |8^ ) the two-rung terms have been subtracted 
in order to avoid a double counting of terms in the self- 
energies. 

We now turn to the analytic continuation to real fre- 
quencies of the expressions derived above. Transforming 



the Matsubara summations into contour integrals shows 
that integrations along branch cuts of auxiliary particle 
Green's functions carry an additional factor exp(— /3A) 
as compared to integrations along branch cuts of phys- 
ical Green's functions, which vanishes upon projection 
onto the the physical Fock space, A ^ oo. Thus, as a 
general rule, only integrations along branch cuts of the 
c-electron propagators contribute to the auxiliary parti- 
cle self-energies. Therefore, by performing the analytic 
continuation, iujn — > w — zO = ix) in all frequency variables, 
we obtain the advanced pseudofermion self-energy. 



(84) 



(c^) = -Af r /-/—/(£- c.) f{e' - uj)Al{uj - e)G,{e) iiN{Q)T!f''\e,e' ,e + e' - u:)Al{u: - e')G,{e') (85) 



and the advanced slave boson self-energy, 
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■^bfi 



(c.) 
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(u;)^-N I -f{e-cu) A<'^{u;-e)nU{0)Tt'\u;,u;,e), 



(87) 



where the vertex functions are given by Eqs. (pO|), (81) and (P2|), ( p3[ ) with 



T 



(c/)(±)p 



t(c.,c.',f^) = c/(^^)w;;(c.,c.',f^) 
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±(+r) / ^f{e^n)GM^ + e-n)G,-,{e)AlAe-^) T^^'^^^ie^io' ,n) 
7r^^i0) U^^'^ ^%iLO, Lo', n) = -r' 1 f f{e - n)Gf,iu; + e - n)G,f,{e)Al,ie - n) Gf^Lu' + e-n) 



(90) 
(91) 



In the above expressions, like in the NCA equations (|5^)- 
(p^), we have used the dimensionless conduction elec- 
tron spectral density, A°{uj) = i ImG°^^(w - iO)/7V(0), 
and we have suppressed obvious spin and channel indices. 
All frequency variables are to be understood as the limit 
oj = uj — iO. 

The equations (|84|)-(p7|), supplemented by the vertex 
fimctions Eqs. (p0|)-(p3^, (p8|)-(^ll) form, together with 



the NCA contributions Eqs. (p2|), ( p3| ) and the definitions 
of the auxiliary particle Green's functions, Eqs. ( [2^ , 
( p3| ) , the closed set of self-consistent CTMA equatioiisEl. 
It is seen that in these equations only those branches 
of the T-matrix vertex functions appear which arc ad- 
vanced with respect to all three frequency variables, al- 
though in general the T-matrix consists of 2^ indepen- 
dent analytical branches. This simplification is a conse- 
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FIG. 7. (1) Pseudofermion and slave boson spectral functions Af and At, in the Kondo regime {N = 2; Ed = —0.05, 
r = 0.01 in units of the half-bandwidth D), for a) the single-channel {M = 1) and b) the multi-channel (M = 2) case. In a) 
the symbols represent the results of NRG for the same parameter set, T — 0. The slopes of the dashed lines indicate the exact 
threshold exponents as derived in section III B for M = 1 and as given by conformal field theory for M — 2. The insets show 
magnified power law regions. (2) CTMA results (symbols with error bars) for the threshold exponents a/ and at of ^/ and 
At, {N = 2, AI = 1). Solid lines: exact values (section III B), dashed lines: NCA results (section V B). 



quence of the exact projection onto the physical sector 
of Fock space. Inspection of the analytically continued 
CTMA equations also shows that the slave boson self- 
energy is obtained from the pseudofermion self-energy, 
including the proper signs, by simply replacing Gj^Gi, 
and inverting the frequency argument of in all expres- 
sions. 

Eqs. (p^)-(^ll) may be rewritten in terms of the spec- 
tral functions without threshold, Aj^f,, in a straight- 
forward way as explained in section V C, thus avoid- 
ing divergent statistical factors in the d-electron Green's 
function. 

The self-consistent solutions are obtained by first 
solving the linear Bethe-Salpeter equations for the T~ 
matrices by matrix inversion, computing the auxiliary 
particle self-energies from T^^f^ and T^^^), and then 
constructing the fermion and boson Green's functions 
from the respective self-energies. This process is iter- 
ated until convergence is reached. We have obtained 
reliable solutions down to temperatures of the order of 



at least 1Q~^Tk both for the single-channel and for the 
two-channel Anderson model. Note that Tk ^ in 
the Kondo limit; in the mixed valence and empty im- 
purity regimes, significantly lower temperatures may be 
reached, compared to the low temperature scale of the 
model. 



C. Results for the auxiliary particle spectral 
functions 

As shown in Fig. Q (1) a), the auxiliary particle spectral 
functions obtained from CTMAe3 are in good agreement 
with the results ot-a numerical renormalization group 
(NRG) calculatioiJ2a (zero temperature results), given 
the uncertainties in the NRG at higher frequencies. Typ- 
ical behavior in the Kondo regime is obtained: a broad 
peak of width ~ T in A;, at w ~ \Ed\, representing the 
hybridizing d-level and a structure va. Af at w ~ Tk- 
Both functions display power law behavior at frequen- 
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cies below Tr-, which at finite T is cut off at the scale 
a; ~ T. The exponents extracted from the frequency 
range T < lo < Tk of our finite T results compare well 
with the exact result also shown (see insets of Fig. 0(1) 
a)). A similar analysis has been performed for a number 
of parameter sets spanning the complete range of d-level 
occupation numbers Ud- The extracted power law expo- 
nents are shown in Fig. |^ (2), together with error bars 
estimated from the finite frequency ranges over which 
the fit was made. The comparatively large error bars 
in the mixed valence regime arise because here spin flip 
and charge fluctuation processes, described by the poles 
in T'^'^^) and T^'^''\ respectively, are of equal importance, 
impeding the convergence of the numerical procedure. 
In this light, the agreement with the exact results (solid 
curves) is very good, the exact value lying within the 
error bars or very close in each case. 
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FIG. 8. (1) Pseudofermion and slave boson spectral func- 
tions Af and At for the X-ray problem (A^ = 1, M = 1; 
Ed = 0.005, r = 0.003 in units of the half-bandwidth D) 
[48] . The slopes of the dashed lines indicate the exact thresh- 
old exponents as derived in section III B. (2) CTMA results 
(symbols) for the threshold exponents a / and at of the spec- 
tral functions Af and At, {N = 1, M = 1). Solid lines: exact 
values (section III B). 

In the multi-jciiannel case {N > 2, M > N) NCA 
has been shownE3 to reproduce asymptotically the cor- 
rect threshold exponents, a/ = M/{M + N), — 
N/{M + N), in the Kondo limit. Calculating the T- 
matrices using NCA Green's functions (as discussed in 



the single-channel case) we find again a pole in the sin- 
glet channel of T^'^^\ However, in this case the CTMA 
does not renormalize the NCA exponents in the Kondo 
limit of the two-channel model, i.e. the threshold expo- 
nents obtained from the CTMA solutions arc very close 
to the exact ones, af ~ 1/2, af, ~ 1/2, as shown in Fig. 
(1) b). 

We note in passing that we have also solved the spin- 
less caS|Q-|(A^ = 1, M = 1) using the auxiliary particle 
methodeS. The Anderson impurity model, Eqs. (||)-(H), 
then reduces to a noninteracting resonant level model. 
In this sense, the spinless case may be seen as the most 
extreme case of a quantum impurity model with a Fermi 
liquid ground state. It does, however, remain non-trivial 
in the slave particle representation Eq. ( pd] ) and, there- 
fore, constitutes a test case for the description of Fermi 
liquid behavior within the auxiliary particle method. As 
discussed in section III B, the auxiliary particle spectral 
functions Af, Af, in general display an infrared thresh- 
old with power law behavior which is induced by an 
orthogonality catastrophy analogous to that of the X- 
ray problem. In the spinless case, Af and At have been 
shown to correspond precisely to the photoemission and 
X-ray absorptioit-spectral densities of the X-ray prob- 
lem, respectivelylla. The results of the CTMA for the 
iV = 1, M = 1 case are shown in Fig. ^ Good quanti- 
tative agreement of the X-ray exponents extracted from 
the numerical evaluation of the CTMA equations is found 

in the regions ~ 1 and ^ which may be identi- 
fied with the regions of a strong (scattering phase shift 
~ 7r/2) and a weak potential (scattering phase shift 
ijcr — 0) scatterer. 

The agreement of the CTMA exponents with their ex- 
act values in the Kondo, mixed valence and empty im- 
purity regimes of the Spin 1/2 single-channel model, in 
the spinless model and in the Kondo regime of the two- 
channel model may be taken as evidence that the T- 
matrix approximation correctly describes both the Fermi 
liquid (iV = 1,2, M — 1) and the non-Fermi liquid 
(N = 2, A/ = 1,2) regimes of the SU(N)xSU(M) An- 
derson model. Therefore, we expect the CTMA to cor- 
rectly describe physically observable quantities of the 
SU(N)xSU(M) Anderson impurity model as well. 



D. Results for physical quantities: Spin 
susceptibility 

We have calculated the static spin susceptibility x of 
the Anderson model in the Kondo regime by solving the 
CTMA equations in a finite magnetic field H coupled to 
the impurity spin and taking the derivative of the mag- 
netization M = ^.9Ms(?i/T~"-/i) with respect to H. The 
resulting x(r) — {dM/dH)T is shown in Fig. |^ both for 
the single-channel case {N — 2, M — 1) and for the 
two-channel Anderson Model (A^ = 2, M = 2). It is 
seen that in the single-channel case CTMA correctly re- 
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produces the constant Pauli susceptibility (Fermi liquid 
behavior) below Tk, while NCA gives an incorrect, non- 
analytic temperature dependence of x(T) — x(0) oc —T'^/'^ 
at low T. In the two-channel case CTMA des<;ribcs the 
non-Fermi liquid behavior, i.e. the the exactll3 logarith- 
mic temperature dependence of the susceptibility below 




iogio(T/TK) 



the Kondo scale Tk- In contrast, the NCA solution re- 
covers the logarithmic behavior only far below Tk- Other 
physical quantities have been calculated for the Anderson 
model at low temperaJtures as well and will be published 
in forthcoming worlJla. 




IoQioCT/Tk) 



FIG. 9. Static susceptibility of the single-channel (A'^ = 2, M = 1) and the two-channel (A'^ = 2, M = 2) Anderson impurity 
model in the Kondo regime {Ed = — O.SD, F = O.ID, Lande factor g = 2). In the single-channel case, the CTMA and NCA 
results are compared with the Bethe ansatz result for the Kondo model [14], where the CTMA as well as the Bethe ansatz 
curves are scaled by their respective Kondo temperatures [14]. The CTMA susceptibility is universal in that it depends only 
on the ratio T/Tk for various values of the parameters of the Anderson model (not shown). 



VII. ANDERSON MODEL AT FINITE U: 
GENERALIZED NCA AND CTMA 

The single channel Anderson model at finite U may be 
represented in terms of the pseudofermion operators fa 
and slave boson operators a, h (discussed in the beginning 
of section III) as 



H = H, + Ed{^af+U + 2a+a) + Ua+a 



ka 



ctJb+U+VafIaa) + h-C- 



subject to the constraint 



(92) 



(93) 



There are now two bosons, a "light" boson b (as in the 
case U ^ oo considered before) and a "heavy" boson 
a. As seen from Eq. (p2[), the strong two-particle in- 
teraction U is transformed into a potential term by this 
representation. The corresponding diagrams are made of 
propagator lines for the light boson (wiggly line, as be- 
fore) and for the heavy boson (zig-zag line) as well as 
a new vertex, where incoming conduction electron lines 
and pseudofermion lines merge into an outgoing heavy 
boson line. 



A. Generating functional 

It seems straightforward to define a generalization of 
NCA for finite U by adding to the second order skele- 
ton diagram for the generating functional $ a second 
one where the light boson is replaced by the heavy boson 
(see Fig. ^ a)). Ste approximation has been consid- 
ered sometime agoEJ'Ell, where it was found to fail badly: 
Not even the Kondo energy scale is recovered in the so 
defined approximation. The reason for this failure is ob- 
vious: In the Kondo regime (n^ ^ 1) the local spin is 
coupled to the conduction electron spin density at the im- 
purity through the antiferromagnetic exchange coupling 



\ ad ; Ed+U' 

this relation arise from virtual transitions into the empty 
and doubly occupied local level, which contribute equally 
in the symmetric case | |= Ed + U. The symmetric 
occurrence of these two virtual processes in all intermedi- 
ate states is not included in the simple extension of NCA 
proposed above. The self-energy insertions in each of the 
two diagrams contain always only one of the processes, 
leading to an effective J which is only one half of the 
correct value. Correspondingly, the Kondo temperature 
Tk ^ exp(— l/A/'(0) J) comes out to scale as the square 
of the correct value, which is orders of magnitude too 
small. 



). The two terms on the r.h.s. of 
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FIG. 10. Diagrammatic representation of the generating functional for the Anderson impurity model at finite U. Solid, 
dashed, wiggly and zig-zag lines correspond to conduction electron c, pseudofermion /, "light" boson b, and "heavy" boson a 
propagators, a) NCA including light (empty impurity) and heavy (doubly occupied impurity) boson lines, b) UNCA containing, 
in addition to an infinite number of (bare) hybridization processes into the empty impurity state (non-skeleton diagrams), a 
single hybridization into the doubly occupied impurity state, c) Symmetrized, finite-?7 NCA (SUNCA): The infinite series of 
all diagrams of a)-c) generates an approximation where the hybridization processes into the empty and doubly occupied state 
are treated in a symmetric way, i.e. for each contribution with a bare light boson line there is a corresponding contribution 
with the bare light boson line replaced by a bare heavy boson line, and vice versa. 



To correct this deficiency it is necessary to include ad- 
ditional diagrams, restoring the symmetry between the 
two virtual processes, excitation into the empty or into 
the doubly occupied impurity state. Due to the exact 
projection onto the physical Hilbert space, each diagram 
contributing to <f> contains exactly one auxiUary parti- 
cle ring (see, e.g. Fig. ||). As a first step one may add 
the next order diagram to $ (see Fig. |l^ b)). As will 
be seen below, this approximation, later referred to as 
UNCA, helps to recover a large part of the correct behav- 
ior of Tk- a completely symmetric treatment, however, 
requires the summation of infinite classes of diagrams. 
These diagrams are generated by replacing a light boson 
line with a heavy boson line in each of the bare (non- 
skeleton) diagrams of the elementary NCA diagrams of 
Fig. |ll|. Each replacement leads to a crossing of conduc- 
tion electron lines spanning one fermion and at most two 
boson lines. As a result, the class of diagrams of <I> shown 
in Fig. ^ c) are obtained. These diagrams look similar 
to the CTMA diagrams shown in Fig. ||, but contain one 
light boson line and an arbitrary number of heavy bo- 
son lines, or vice versa. Diagrams with, for example, two 
light boson lines and an arbitrary number of heavy boson 
lines (and conduction electron lines spanning at most one 
fermion line) are reducible and do not appear. We will 
call the approximation defined by the generatin g fu nc- 
tional given by the sum of the diagrams of Figs. ^ a)- 
c) "symnrfitpzed, finite-C/ non-crossing approximation" 
(SUNCA)l3e3. It is interesting to note that the sym- 



metric treatment of light and heavy bosons leads to dia- 
grams for $ of a structure similar to the ones defining the 
CTMA (Fig. H). The above approximation corresponding 
to the CTMA at J7 ^ oo, termed "symmetrized, finite-C/ 
Conserving T-matrix Approximation (SUCTMA) is thus 
defined in a natural way by summing up all skeleton dia- 
grams with an arbitrary number of light or heavy boson 
lines, dressed by conduction electron lines spanning only 
one fermion line. The rational for keeping only these 
crossings is that (as explained in the appendix). 




FIG. 11. Two examples of how diagrams involving hy- 
bridization into the doubly occupied impurity state are gen- 
erated from the bare noncrossing diagrams of the infinite U 
case by replacing light with heavy boson lines (see text). 
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diagrams containing conduction electron lines spanning 
two or more fermion lines may be grouped into sets 
where the most singular contributions cancel. Thus, the 
SUCTMA is defined by adding to the diagrams of the 



SUNCA the CTMA diagrams (Fig. |) with (only) light 
boson lines or (only) heavy boson lines. The SUCTMA 
equations have not yet been evaluated. 



a) 




+ 




b) 




+ 




FIG. 12. Diagrammatic representation of the Bethe-Salpeter equations for a) the renormalized light boson (empty impurity) 
and b) the renormalized heavy boson (doubly occupied impurity) vertex, as generated by the SUNCA Luttinger-Ward functional 
(Fig. 11). 



B. Results of SUNCA 

The pseudoparticle self-energies , and E/ as well 
as the local conduction electron self-energy Sc are ob- 
tained as funcional derivatives with respect to the cor- 
responding Green's functions. They can be expressed in 
terms of the two vertex functions defined in Fig. |l^ as 
solutions of Bethe-Salpeter equations. It is seen that the 
vertex functions are obtained by ladder summations in- 
volving a heavy boson or a light boson line as the rung 
of the ladder. These three-point vertex functions depend 
on two frequencies and are much less difficult to calcu- 
late than the four-point T-matrices of the CTMA. The 
full expressions of the self-energies are given in Fig.[l^. 
The local d-electron Green's function after projection is 
again proportional to T,c'. 



U = 



£f = 



V^Gh = 



1 

U2 



(94) 



FIG. 13. Diagrammatic representation of the auxiliary 
particle selfenergies of SUNCA in terms of the renormalized 
hybridization vertices, defined in Fig. 12. Note that in each 
line the third diagram is subtracted in order to avoid double 
counting of terms within the first two diagrams. 



The self-consistent set of equations defined by Fig. |l3| 
and the expression for the pseudo-particle Green's func- 
tion, Eq. (pl|), supplemented by the expression for the 
heavy boson Green's function 



Qa{it^n) = \i^n " (2£^d + U) - Ea(iw„)| 



(95) 



have been solved numericallyEEEi. As expected, the pseu- 
doparticle spectral functions are found to display power 
law divergencies at the infrared threshold. The exact 
power law exponents, derived from the Friedel sum rule 
argument based on the assumed Fermi liquid ground 
state (section III B), are given by 



"3 



-1 



(96) 
(97) 
(98) 



The threshold exponents obtained from the SUNCA so- 
lution are shown in Fig. |lj as functions of the aver- 
age impurity occupation number rid. They agree with 
the exact results in the Kondo limit {ud 1) where 
aj = ab = aa = ^, but deviate for < 1. It is ex- 
pected that the SUCTMA will recover the correct expo- 
nents, as does the CTMA in the case of infinite U. The 
d-electron spectral function for the symmetric Anderson 
model as obtained in SUNCA is shown in Fig. |l^ (left 
panel). In this case the Kondo resonance is located ex- 
actly at lli — 0. Its width is a measure of the Kondo 
temperature Tr-. The three curves shown correspond to 
the elementary NCA (Fig. [lO|a)), the UNCA (Fig. |l^ b)) 
and the SUNCA. The inset shows that for the parameters 
chosen the Tk obtained in NCA is too small by a 
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FIG. 14. Infrared threshold exponents of the auxiliary par- 
ticle spectral functions for the case of finite U in dependence 
of the the impurity occupation Ud- Dashed curved lines: ex- 
act results (Eqs. (96)-(98)); horizontal lines: NCA results; 
data points with error bars: SUNCA results. In the Kondo 
limit {ud — » 1) the exact exponents are recovered, while in 
the mixed valence and empty impurity regime the SUNCA 
results for af and at cross over to the NCA values. 

factor 10~^, whereas-the UNCA is only a factor of 1/3 
off. Contrary to Ref.E3 we do not find that the UNCA is 
sufficient to recover Tk- In Fig. ^ (right panel) the re- 
sults of the SUNCA for Tk as determined from the width 
of the Kondo resonance are compared to the exact result, 



Tk — min 



(99) 



where / — 



2ru 



. The agreement is seen to be very 



good for a large range of parameters / and Ed- 



VIII. CONCLUSION 

We have reviewed a technique to describe correlated 
quantum impurity systems with strong onsite repulsion, 
which is based on a conserving formulation of the aux- 
iliary boson method. The conserving scheme allows to 
implement the conservation of the local charge Q with- 
out taking into account time dependent fluctuations of 
the gauge field A, while exactly projecting the quantum 
dynamics onto the physical subspace of no double oc- 
cupancy of sites (in the limit of infinitely strong on-site 
repulsion). Any spurious condensation of the auxiliary 
boson field is avoided in this way. 

By including the physically dominant contributions, 
spin flip processes in the Kondo regime, and spin as 
well as charge fluctuation processes in the mixed va- 
lence and empty impurity regimes, the method recovers 
the Fermi liquid ground state of the single-channel An- 
derson impurity model as well as the non-Fermi-liquid 
low-temperature behavior of the two-channel Anderson 
model: The correct infrared threshold exponents of the 
auxiliary particle propagators, which are identified as in- 
dicators for Fermi or non-Fermi liquid behavior, are ob- 
tained in both cases. Physical quantities, like the mag- 
netic susceptibility, are correctly described both in the 
Fermi and in the non-Fermi liquid cases of the model over 
the complete temperature range, including the crossover 
to the correlated many-body state at the lowest temper- 
atures, which has previously been notoriously difficult to 
obtain within a diagrammatic many-body theory. 




FIG. 15. Left panel: Local electron spectral function calculated using NCA, UNCA, and SUNCA. The Kondo temperature 
is determined by the peak width. It is seen that in NCA it comes out orders of magnitude too low. — Right panel: Kondo 
Temperature for various parameters Ed, U and fixed F. Solid lines represent the results from the Bethe ansatz solutions, Eq. 
(99). Data points are the SUNCA results determined from the width of the Kondo peak in the d-electron spectral function. 
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We have also presented a generalization of the auxiliary 
particle method for finite on-site repulsion U and solved 
the corresponding equations on the level of a generalized 
non-crossing approximation (SUNCA) , where virtual ex- 
citations into the empty or into the doubly occupied im- 
purity state are treated in a fully symmetric way. In this 
way the correct low-temperature scale of the finite-U An- 
derson model was obtained. The controlled treatment of 
a finite on-site repulsion is a precondition for applying 
the auxiliary particle method to models of the Hubbard 
type or - within the scheme of the dynamical mean field 
theory (limit of infinite dimensions) - of the Anderson 
impurity type, where a Mott-Hubbard metal-insulator 



transition occurs as a function of band filling for finite 
[/. As a standard diagram technique the conserving aux- 
iliary particle method has the potential to be applicable 
to problems of correlated systems on a lattice as well as 
to mesoscopic systems out of equilibrium via the Keldysh 
technique. These developments are currently under in- 
vestigation. 
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APPENDIX: INFRARED CANCELLATION OF 
NON-CTMA DIAGRAMS 

The CTMA is not only justified on physical grounds 
by the inclusion of the maximum number of spin flip and 
charge fluctuation processes at any given order of pertur- 
bation theory, but also by an infrared cancellation of all 
diagrams not included in the CTMA. In the following we 
will prove this cancellation theorem. 



(1) Power counting. — Each auxiliary particle loop car- 
ries a factor of the fugacity exp(— /3A), which vanishes 
upon projection onto the Q = 1 subspace, A ^ oo. 
Therefore, an arbitrary / or 6 self-energy diagram con- 
sists of one single line of alternating fermion and boson 
propagators, with the hybridization vertices connected 
by conduction election lines in any possible way, as shown 
in Fig. |l^ (see alsocS). Such a fermion self-energy skele- 
ton diagram of loop order L is calculated as 
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where Gf^b are the renormalized , i.e. power law divergent 
auxiliary particle propagators, and Lgp and Lch denote 
the number of spin (or fermion, c-f) loops and the num- 
ber of channel (or c-b) loops contained in the diagram, 
respectively. Spin and channel indices are not shown 
for simplicity. Each of the auxiliary particle frequencies 
LUi, Lo[ coincides with one of the integration variables e^, 
j = 1, . . . , L, in such a way that energy is conserved at 
each hybridization vertex. This implies that the sign of 
the frequency carried by a c-electron line is Si = +, if 
the c-electron line runs from right to left, and = — , 
if it runs from left to right in Fig. |l^. An analogous ex- 
pression holds for the slave boson self-energy diagrams. 




FIG. 16. Typical pseudofermion self-energy skeleton di- 



agram of loop order L 
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By substituting Xj — Sj/uj, j — 1, . . . ,L and factoring 
out oj^'^f-C') from each fermion (boson) propagator, the 
infrared behavior of the term Eq. (|A1|) is deduced as 



(A2) 



where C is a finite constant. Clearly, when the NCA solu- 
tions are inserted for the propagators G/,b, i.e. af + at = 
1, their power law be havi or is just reproduced by any 
term of the form Eq. (Al). However, this is no longer 
the case for the exact propagators in the Fermi liquid 
regime (M < N), where in general af + ah > 1. Thus, 
the infinite resummation of terms to arbitrary loop order 
is unavoidable in this case. 

(2) Infrared cancellation. — As discussed in section VI 
B, the CTMA is equivalent to the self-consistent summa- 
tion of all skeleton free energy diagrams, where a conduc- 
tion electron line spans at most two hybridization ver- 
tices (Fig. 1^). Thus, any skeleton self-energy diagram 
not included in CTMA contains at least one conduction 
electron "arch" which spans four (or more) vertices, with 
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FIG. 17. Set of contributions to skeleton diagrams not 
contained in CTMA which cancel in the infrared limit to lead- 
ing and subleading order in the external frequency, — > 0. 



four conduction lines reaching from inside to outside of 
the arch as shown in Fig. 0(a). For each such dia- 
gram there exists another skeleton, which differs from 
Fig. ^(a) only in that the end points of two conduc- 
tion hues inside the arch arc interchanged (Fig. [l^(b)). 
The corresponding permutation of fermionic operators 
implies a relative sign between the terms Fig. ^(a) and 
(b) . Without loss of generality we now assume w > for 
the external frequency of the self-ener gy. The leading 
infrared singular behavior of the term (Al) arises from 
those parts of the integrations, where the arguments of 
the Gf,Gb are such that the divergences of all propaga- 
tors lie within the integration range. This implies at least 
—UJ < Sj < 0, j = 1, . . . , L. Therefore, the terms corre- 
sponding to Fig. 0(a), (b) differ only in the frequency 
arguments of the Green's functions inside the arch, and 
at temperature T = the leading infrared behavior of 
their sum reads. 
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F{uj, {sj}) [Gf{ijj + Lo' + uji)Gb{ijj + uj' + UJI+ UJ2) - Gf{Lo + uj' + Lo-i)Gb{uj -I- cj' + -f ^2)] 



Here ui denotes the sum of all frequencies to, ej enter- 
ing the diagrammatic part. Fig. from the left, and 
F{llj, {ej}) consists of all terms which are not altered by 
interchanging the c-electron lines. In the infrared limit, 
UJ1—UJ3 0, the term in square brackets may be written 
as 

[Gf{uj + Lu')Gb{uj + Lu' + UJ2)] (wi - LU3) (A4) 

and upon performing the integrations over lui, W3 the 
difference (wi — 013) leads to an additional factor of uj. 
A similar cancellation of the leading infrared singularity 
occurs between the terms shown in Fig. |l^(c), (d). In 
an analogous way it may be shown that combining the 
terms Fig. ^ (a)-(d) leads to a factor of compared 



I 

to the power counting result for one single term. Thus, 
the infrared singularity of all non-CTMA terms of loop 
order L is weaker than the Lth order CTMA terms by at 
least 0(^2), 

+ ■ • • + ^f'''\Lu) ^"/..+i(l-"/-".)+2. 

(A5) 

It should be emphasized that in the above derivation, L 
appears only as a parameter and, thus, the cancellation 
theorem holds for arbitrarily high loop order L. This 
proves that the CTMA captures the leading and sub- 
leading infrared singularities (uj 0) at any given order 
L. 
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